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The one-step replica symmetry breaking (RSB) is used to study a two-sublattice fermionic infinite- 
range Ising spin glass (SG) model in a transverse field F. The problem is formulated in a Grassmann 
path integral formalism within the static approximation. In this model, a parallel magnetic field 
H breaks the symmetry of the sublattices. It destroys the antiferromagnetic (AF) order, but it 
can favor the nonergodic mixed phase (SG-I-AF) characterizing an asymmetric RSB region. In this 
region, intra-sublattice disordered interactions V increase the difference between the RSB solutions 
of each sublattice. The freezing temperature shows a higher increase with H when V enhances. A 
discontinue phase transition from the replica symmetry (RS) solution to the RSB solution can appear 
with the presence of an intra-sublattice ferromagnetic average coupling. The F field introduces a 
quantum spin flip mechanism that suppresses the magnetic orders leading them to quantum critical 
points. Results suggest that the quantum effects are not able to restore the RS solution. However, 
in the asymmetric RSB region, F can produce a stable RS solution at any finite temperature for 
a particular sublattice while the other sublattice still presents RSB solution for the special case in 
which only the intra-sublattice spins couple with disordered interactions. 

PACS numbers: 75.10.Nr, 75.10.Jm, 64.60.Cn 



I. INTRODUCTION 

The quantum spin glass (SG) theory^ has been sub- 
ject of extensive studies mainly by its conceptual diffi- 
culties to describe the ordered phase. This problem be- 
comes even more interesting when the SG and a long 
range antiferromagnetic (AF) interaction are both con- 
sidered. The glassy AF order can be source of complex 
physics processes. For instance, quantum fluctuations 
and the interplay among disorder, frustrated magnetism 
and other interactions have been an essential feature 
of some heavy fermion systems^ and high temperature 
superconductorsi^ It can lead to a non-Fermi liquid be- 
havior and to a quantum critical point (QCP). In this 
case, it is important to establish a theoretical framework 
to treat the SG/AF competition with other physical pro- 
cesses such as Kondo effect^ or superconductivitjs^ in the 
same footing at low temperature in which quantum ef- 
fects become significant. We do not intend discussing 
this specific problem, but we concentrate on the quan- 
tum magnetic properties of a fermionic formulation for a 
two-sublattice transverse field Ising SG model. 

The infinite range Sherrington-Kirkpatrick SK model^ 
in a transverse tunneling field F has been extensively 
studied as a quantum SG model that presents a QCP and 
an experimental realization in the system LiIl0a;Yi_3;F4ii 
However, even this simplest theoretical quantum SG 
model has recently been theme of controversial results 
concerning the local stability of replica symmetry (RS) 
solution (see Refs. 8,9,10,11,). For instance, the classical 
SK model is well described by Parisi's replica symmetry 
breaking (RSB) solution. It considers the SG phase as a 
nonergodic situation in which the free energy landscape is 
characterized with many valleys separated by large barri- 
ers. This raises the question of whether or not the quan- 



tum tunneling through these barriers is able to restore 
the RS when quantum effects are included (F > 0). 

Concerning the SG/AF competition, a two-sublattice 
Ising spin model with SK-type interactions between clas- 
sical spin variables of distinct sublattices has been in- 
troduced by Korenblit and Shender (KS) to discuss the 
SG/AF problem at mean field level>i^ A generalization 
of the KS model (GKS model) allows inter- and intra- 
sublattice disordered interactions.^'''-- These two models 
(KS and GKS) have shown an interesting feature related 
to a parallel magnetic field H. It breaks the symmetry 
between the sublattices that can favor frustration increas- 
ing the RSB temperature rangej^s^ii^ The RSB with 
this asymmetry caused by the H field is named here as 
asymmetric RSB region, in which the SG and the usual 
AF order coexist (mixed phase). However, a quantum 
version for the GKS model has received less considera- 
tion. In this case, the quantum effects on the asymmetric 
RSB region remain an open question. For instance, the 
RSB occurs at the same temperature in both sublattices 
for the classical GKS modeli^ However, the degree of 
RSB can be distinct in each sublattice within the asym- 
metric RSB region. Regarding that difference, we can 
ask how the RSB behavior is when quantum fluctuations 
are present in the GKS model. 

On the other hand, there is a particular situation for 
the classical GKS model (PGKS) in which the inter- 
sublattice interactions are only antiferromagnetic with- 
out disorder and the intra-sublattice spins only couple 
with disordered interactions.^^ In this case, the RSB can 
first occur for a specific sublattice (let us say b), while 
the other sublattice (a) still presents RS solution. The 
sublattice a can also show a RSB solution, but at a lower 
temperature. Now, the quantum tunneling is more pro- 
nounced at lower temperature. In this context, the RSB 
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solution of the sublattice a can be more affected by the 
quantum fluctuations than the sublattice b. In this case, 
could the sublattice a present a stable RS solution for 
any finite temperature at the same time that the sublat- 
tice b still shows a RSB solution when F > 0? That is 
an important point for the SG/AF problem because the 
solution of a particular sublattice also depends on the 
order parameters of the other sublattice.^-:-'* Moreover, 
the frustration in the asymmetric RSB region could be 
changed if only the spins of the sublattice b present RSB 
solution. 

The aim of this paper is to investigate the effects of 
a transverse magnetic field F on a fermionic formulation 
for the GKS model where inter- and intra-sublatticc SK- 
type interactions are considered. In special, we investi- 
gate the effects produced by the F field on the frustration 
in the asymmetric RSB region for the case in which only 
the intra-sublattice interactions are disordered (PGSK 
model). The Parisi's scheme^^ of one-step RSB (IS-RSB) 
is used to study the non-trivial ergodicity breaking re- 
gion. In this case, we check if the quantum fluctuations 
are strong enough to cause tunneling between the free 
energy barriers of the many degenerated thermodynamic 
states. 

The H field breaks the symmetry between the sub- 
lattices in the GKS modelJ^ As a consequence, H can 
favor frustration at the same time that it destroys the 
AF order, so the freezing temperature Tf can increase 
with H. On the other hand, the F field introduces a 
quantum spin flipping mechanism, which can suppress 
the magnetic orders leading their critical temperatures to 
QCPs.^^ Therefore, both magnetic fields produce effects 
that can compete on Tf in this two-sublattice SG/AF 
problem. For example, in a previous fermionic version 
for the KS model in a transverse field (without intra- 
sublattice interaction),— results within the RS solution 
have shown that the effects produced by F are dominant 
at low H. However, when H enhances, the asymmetry 
between the sublattices becomes important and Tf in- 
creases. This scenario can be even more complex when 
the intra-sublattice disordered interactions are consid- 
ered. They could affect the Tf behavior in the asymmet- 
ric RSB region. Furthermore, a first-order phase tran- 
sition from the AF to the paramagnetic phase can raise 
when H increasesii^ Therefore, in this paper, we also 
study the interplay between disorder and magnetic fields 
{H, F) in the SG/AF competition for the fermionic for- 
mulation of the GKS model in a transverse field within 
IS-RSB solution. 

In this work, the hierarchical tree-structure is adopted 
in each sublattice to obtain a RSB solutionJ^ Thus, the 
elements of the replica matrices for the sublattices a and b 
are parameterized following the Parisi scheme of IS-RSB. 
These replica matrices are independent of each other for 
the PGKS model in which the two sublattices present 
independent IS-RSB solutions. However, the replica ma- 
trices for the sublattices a and b are coupled for the 
GKS model. In this case, we adopt the approximation 



nia = rrib in order to decouple them, where vrip {p = a ot 
b) is the Parisi block-size parameter-i? for the sublattice 
p. In this approximation, the situation nia 7^ m-b never 
occurs. Nevertheless, it still allows sublattice dependence 
for the behavior of the other IS-RSB order parameters. 
An alternative RSB method based on the Parisi's scheme 
and a wave-like structure has been proposed in Ref. [l^ 
for the KS model without intra-sublattice disorder. It 
is the modulated RSB scheme that considers a possible 
asymmetry between rua and mf, which is controlled by 
a variational parameter in the replica limit. However, 
this RSB scheme increases the numerical efforts to ob- 
tain the IS-RSB results and it can only provide a very 
small improvement in the free energy compared with the 
Parisi scheme. Therefore, we do not use the modulated 
RSB scheme in the present GSK model, in which the 
intra-sublattice SK-type interactions add some difhcul- 
ties when this RSB scheme is implemented. 

The partition function of the present work is computed 
with the path integral formalism in which the spin oper- 
ators are represented by bilinear combinations of Grass- 
mann variables as in Ref. J^7. The IS-RSB free energy is 
calculated within the static approximation (SA) that is 
used to treat the spin-spin correlation functions.* The SA 
deserves some comments. It is reasonable to obtain phase 
boundaries even at low temperature.— However, it can 
yield inaccurate quantitative results for the dynamical 
behavior of the correlation functions at very low temper- 
ature when F > 0. Nevertheless, there are several recent 
paper a*'^'**!*^'*^ adopting this particular approximation 
because of a lack of a suitable analytical method to go 
beyond the SA in the whole SG phase in the femionic for- 
malism, in which the SG phase has also a shortcoming 
related to the RSB solution. Furthermore, when temper- 
ature decreases for small and intermediary values of F in 
the one-lattice SG problem,** the results produced by SA 
with IS-RSB are in good qualitative agreement with that 
obtained by Trotter-Suzuki formalism in which numerical 
methods are used to treat the dynamic of the spin-spin 
correlations.^ Moreover, a Landau-Ginzburg formalism 
has been proposed in Ref. [13 to analyze the SG of quan- 
tum rotors near the QCP. In this formalism, the RSB 
occurs in the SG phase at finite temperature T and it is 
suppressed when T — > for the region close to the QCP. 
The results with SA agree qualitatively with this RSB 
behavior when F increases towards the QCP, in which 
the RSB effects become less pronouncedJJ- 

This paper is structured in the following way. In the 
next section, the model is introduced, and the IS-RSB 
free energy is obtained. In section Hill numerical results 
for the IS-RSB order parameters, free energy and phase 
diagrams showing the SG/AF competition in tempera- 
ture versus parallel magnetic field are reported for sev- 
eral values of disorder and transverse magnetic field. The 
last section is left for summary and conclusions. 
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II. GENERAL FORMULATION 

Let us consider a set of localized Ising spins at sites 
of two identical sublattices a and b. The spins of same 
and distinct sublattices can interact with infinite-range 
random couplings. Transverse (F) and parallel (H) mag- 
netic fields are coupled with spins of both sublattices. 
The model is represented by the Hamiltonian 



Op 



p=a,b ipjp 



ia jb 



(1) 



where the sums run over the N sites of each sublattice 
p {= a 01 h). The inter- and intra-sublattice exchange 
interactions Ji^j,, and Vi j are random variables that 
follow Gaussian distributions with mean values —AJq/N 
and 2Vo/N and variances ,P/32N and V'^/16N, respec- 
tively. The spin operators in Eq. ([1]) are defined in terms 
of fermion operators^^: 



and 5,^ = ^(4,TC^pi+<iC^pT) (2) 



where fii^cr — c] ^Ci^a gives the number of fermions at 
site ip with spin projection a =] ov [. c\^^ and Ci^cr are 
the fermions creation and annihilation operators, respec- 
tively. In this representation, the operator 5^ has two 
magnetic eigenvalues and two nonmagnetic eigenvalues 
(empty and double occupied site). Therefore, there are 
two states that do not belong to the usual spin space. 
The present fermionic SG/AF problem admits these four 
states. However, the chemical potential /i is chosen to 
ensure the average occupation of one fermion per site.— 
Lagrangian path integral formalism with spin oper- 
ators represented as anticommuting Grassmann field 
(0, (j)*) is used to treat this fermionic problem.— 
Therefore, the partition function is written as: 



Di 



with the action 



p=a,bjp,<J 



^0 



(3) 



(4) 



f3 = 1/T (T is the temperature) and /i = 0, which corre- 
sponds to the half-filling situation. 

In Eq. the Fourier decomposition of the time- 

dependent quantities is employed. The average over 
the quenched disorders is performed by means of replica 
method {I3F ^ -ikj lim((Z")j,. - l)/n) with the stan- 

dard Hubbard- Stratonovich transformations used to de- 
couple different sites. Within the static ansatz, this pro- 
cedure results in an effective single-site problem that in- 
volves saddle-point conditions for each sublattice mag- 
netization Afp, replica diagonal spin-spin correlation r^. 



replica off-diagonal elements q^''^ which are related to 
the SG order parameters, and also fields Q'^ {= {q"^ + 
q'^'')/2) that couple the sublattices. Therefore, the repli- 
cated partition function is 



z(,^) = /?E 

+/3^E 



-JoACA4" + ^E(^^p) 



(5) 



Eine, 



where a and 7 are the replica indices running from 1 to 



%^ j D[<t^'^*rp]e^j>[H;^f] 



(6) 



with 



Hp^^ = E [^M,p + 2/? (VoM; - JoM;,) S'^] 



OCX Q'y 

Op Op , 



(7) 



and p = a {p' = b) or p = b (p' = a). In Eq. ([7|), 



(9) 



^p^w) = (^4>p^{lu) 4)p*{uj)j is a Grassmann spinor, of 

and of are the Pauli matrices, and lu denotes fermionic 
Matsubara frequencies. 

The Parisi's scheme of one-step replica symmetry 
breaking (IS-RSB) is used to parametrize q^'' asi^ 



rp ii a — J 

qpo if I{a/mp) ^ {l/rup) 



(10) 



where I{x) gives the smallest integer which is greater 
than or equal to x. 

For a general case, when J and V are non-zero (GKS 
model), the Parisi block-size parameter nip is assumed 
to be equal to m for both sublattices. This approxima- 
tion is used in order to follow the standard Parisi scheme 
of RSB. However, the situation with J = and V ^ Q 
(PGKS model) is analyzed by taken vria and as inde- 
pendent parameters {vria ^ mi,) in the next section. 

The parametrization (jlOp with m = nia = rni, is used 
to sum over the replica index. It produces quadratic 
terms in Eq. ([7]) that are linearized introducing new 
auxiliary fields in Eq. The functional Grassmann 

integral is now performed, and the sum over the Mat- 
subara frequencies can be evaluated like Ref. 17. This 
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procedure results in the following expression to the IS- 
RSB free energy: 

2F = -JoMaMb + PJ"^ [m{qaiqbi - qaoqbo) + ran 

-qaiqbi] + { Y^'^P + ^ [nT-illi - 9po) 

p—a,b 



+rl - q^i] 



1 



In 4 



(11) 



with 



/oo 
D^pC0sh{(3^) (12) 
-oo 



where Ap = hp + T^, hp ^ hp + axpVp + a2p, hp = H 
VoMp-JoMp' (p^p), 



aip = j2[V^{qpx - qpo) + J^qp', - g^o)], (13) 



o-2p = ^2[V^{rp - qpi) + .P{r^, - g^'JlCp 
+ ^2{V^qpQ + J^qp,a)zp, 



(14) 



and Dx ~ dx e^^ /^/\/27r (a; = Zp, Vp or ^p). The param- 
eters Qpo, <Zpi, Mp, Tp, and m are given by the extreme 
condition of free energy pTjl . 

The IS-RSB solution of free energy (jlip recovers the 
RS one ((jpi = gpo and m = 0) at temperature T above 
the freezing temperature Tf. In this case, the stability 
against transversal fluctuations of the RS solution follows 
from the Almeida- Thouless analysis;^ which is extended 
to the two-sublattice problem— The Tf is given by the 
higher temperature that satisfies the expression: 



where 



[T] ~ V^da){T] - V'^db) - jUadb = (15) 



jD^p-^smhiP^pY 




/ D^p[^ cosHf3y^p) + ^ siiMP^p)] 



Or 



(16) 



with (Tip ~ in hp. 



III. RESULTS 

Numerical calculations for the set of saddle-point equa- 
tions show different types of solutions. These depend 




FIG. 1: Order parameters versus T for H — and two values 
of T: r = (full lines) and V = 0.8 (dashed lines). 



on the values of various parameters of the model. For 
instance, in the absence of parallel magnetic field H, 
the paramagnetic (PM: qpi = qpo = 0, Mp = 0), the 
spin glass (SG: qpi ^ qpQ > 0, Mp = 0) and the an- 
tiferromagnetic (AF: qpi = qpQ > 0, Ma — —Mh > 0) 
solutions can be found. A mixed phase (SG-I-AF) can 
also be obtained when Ma — —Mi, > with RSB 
(?pi 7^ Ipo) solution, which is described by the order pa- 
rameter Sp = qpi — qpo > 0. The presence of H always 
induces the parameters ^pi and qpQ. In this context, the 
PM solution is characterized by qpi — qpo and Ma = Mb. 
The AF solution is found when Ma ^ M^ with qpi — qpQ. 
Particularly, the AF order is given by the staggered mag- 
netization Ms = {Ma - Mb)/2 > and Sp = 0. The SG 
phase occurs when Mg — with RSB solution. The 
mixed phase is found when Mg > with 6p > 0. 

We first present results for H = 0. In this case, the or- 
der parameters as a function of T arc exhibited in Fig. ([1]) 
for a given degree of frustration rj = VV^ + J^/iVo + Jo)- 
In this work, the quantities T, H, F, Vq and Jq are in 
units of y/V^+lP. For numerical purposes \/V^+J^ = 
1 and 7] = 1/1.9. The order parameter Mg appears con- 
tinuously at the Neel temperature when T decreases. 
The RSB occurs at Tf, in which 6p becomes greater than 
zero. At this point, the staggered magnetization is non 
zero, and this (with dp > 0) characterizes the mixed 
phase. The magnetization of each sublattice becomes 
weakly dependent on T below Tf. Furthermore, when 
H ~ 0, the two sublattices exhibit symmetric results for 
the order parameters, which are identical for any set of 
Vq, Jq, V and J that satisfies the same value of ?/. The 
transverse field F decreases the values of the order pa- 
rameters as well as the critical temperatures (see dashed 
lines of Fig. ([T])). Its effects are more pronounced at low 
temperature. However, it is not able to restore the RS 
solution for H = that is in agreement with Refs. and 
[Til . The free energy of IS-RSB solution can be compared 
with the one of the RS solution in Fig. The dif- 

ference between these two approaches increases when T 



FIG. 2: Free energy versus T for H 
insert shows the resuh for F — 0.8. 



and F = 0. The 



decreases from Tf, and the entropy (5 = —dF/dT) of the 
IS-RSB solution is positive for the range of temperature 
analyzed. 

The results for H > can depend on the set of param- 
eters Vq, Jq, V and J. However, it is not our intention 
to do an exhaustive study of all possible configuration 
of interactions for the model ([1]). Nevertheless, we have 
chosen a specific degree of frustration in which the role 
of each parameter is analyzed. In this case, the ratio 
Tn{H = 0)/Tf{H = 0) remains the same for all results. 
For example, we analyze the behavior of the RSB solu- 
tion and the phase diagrams when the intra-sublattice 
ferromagnetic average coupling Vq is changed, or when 
the intra-sublattice disorder {V) increases. Moreover, the 
main focus of the work is to study the effects of quan- 
tum fluctuations on the asymmetric RSB region. For this 
purpose, we consider two types of coupling: the first one 
(GKS model) admits SK-type interactions between spins 
of same and distinct sublattices (J, V 0). The second 
one (PGKS model) considers only an antiferromagnetic 
coupling (J = 0, Jq > 0) between the sublattices with 
intra-sublattice SK-type interactions {V ^ 0). 



A. Asymmetric RSB for V" / and J / (GKS) 

Results of the IS-RSB order parameters are discussed 
for the first type of coupling in the asymmetric RSB re- 
gion. To do so, we study three distinct configurations 
of inter- and intra-sublattice interactions that satisfy the 
same ry: (i) J ^ V and Vq = Jq = 0.95, (ii) J = V, 
Vo = 1.56 and Jq = 0.34, and (iii) J = 0.24, V = 0.97 
and Vq = Jq — 0.95. In these three cases, we can analyze 
the main effect cased by increasing Vq/ Jq (consider cases 
(i) and (ii)) and by increasing V/ J (consider cases (i) and 
(iii)). 

The i?-dependencc of the order parameters gpi, Spi, 
Mp and Ms is reported in Fig. 1^ for T = 0.3 and F = 0. 
The sequence of panels (a), (b) and (c) in Fig. ^ rep- 
resents the cases (i), (ii) and (iii), respectively. The H 
field breaks the symmetry between the sublattices in the 
region where Mg > 0. However, the results for both sub- 
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FIG. 3: Order parameters versus H for T = 0.3 and F = 0. 
The results are for the following parameters: (a) J = V and 
Vq = Jo = 0.95, (b) J = V,Vo = 1.56 and Jo = 0.34, and (c) 
J = 0.24, V = 0.97 and Vo = Jo = 0.95. 



lattices are symmetric at high values of H {Mg = 0). An 
important point is the i/-dependence of Sp. For instance, 
Sb can increase with H in the mixed phase, which char- 
acterizes an enhancing of the RSB. This Si, odd behavior 
can be associated qualitatively with the 7?-dependence 
of the average internal field hi, (= H + VqMi, — JqMo) 
that acts on spins of sublattice h. For example, there is 
a region at small H where the absolute value of magne- 
tization |Affc| decreases with 77, while Ma has a weaker 
dependence on H than |Mfc| (see Fig. ([3])). As a conse- 
quence, I /if, I decreases with H and 6i, enhances. Now, hp 
is also related to Vq/Jq- When Vq increases as in case 
(ii) (see Fig. ([Sl-b)), the region with solution > de- 
creases. The order parameters can show a discontinuity 




FIG. 4; Order parameters as a function of T for H — 0.5 and 
three values of F (=0.0, 0.8 and 1.0). The interactions are 
given by case (i); J = V and Vo = Jo = 0.95. 



when H crosses from the Mg > to the Mg — solu- 
tion, where dp jumps to a higher value. On the other 
hand, the relations in hp are not directly affected by the 
intra-sublattice disorder V. However, 6}, is higher within 
the mixed phase when V increases as it is shown by com- 
paring cases (i) and (iii) in Figs. ([3]-a) and ^c). For 
instance, the random internal field hp, which acts on the 
sublattice p, also depends on the variance aip (see Eq. 
(fT5)) ). which enhances (decreases) for the sublattice b (a) 
in the asymmetric RSB region when V/ J increases. Con- 
sequently, 5b (da) can increase (decrease) with V/J (Fig. 
(ISl-c)). Therefore, the behavior of the IS-RSB solution is 
also affected by the relation V/ J. In this case, the differ- 
ence between the degree of RSB of the two sublattices is 
increased with V/J, but the RSB always occurs in both 
sublattices simultaneous for J ^ 0. In the next section, 
this result can be compared with the limit case J = 0. 




FIG. 5: Order parameters versus T for H = 0.4, F = 0.0 
and Jo = Vo = 0.95, J = and 1/ = 1.0. The dashed lines 
represent results for the IS-RSB solution when m = ma ~ nib. 
The insert shows the free energy versus T. 



Figure Q shows effects of quantum tunneling on the 
asymmetric RSB region. It exhibits the order parame- 
ters versus T for the case (i) with H = 0.4, F = 0.0, 0.8 
and 1.0. The F field decreases the values of order pa- 
rameters 6a and Sb, but the RSB still occurs at the same 
critical temperature for both sublattices. It means that 
the quantum fluctuations act with the same intensity in 
the RSB solution of sublattices a and b, which still re- 
main asymmetric. Furthermore, the results suggest that 
the quantum tunneling is not able to restore the RS solu- 
tion of neither of the sublattices, even in the asymmetric 
RSB region. The quantum effects on the order parame- 
ters for cases (ii) and (iii) are qualitatively the same as 
those obtained in case (i). 



B. Asymmetric RSB for V" / and J = (PGKS) 

At this point, we restrict our study to the second type 
of coupling, in which the inter-sublattice exchange in- 
teractions are only antiferromagnetic without disorder 
and the intra-sublattice interactions are still disordered 
(Jo ^ 0, J = and V ^ 0). In this case, the structures 
of the replica matrices {Qa} and {Qb} are independent 
of each other. The Parisi's scheme of RSB is applied in 
the two sublattices without any approximation on mp. 
Therefore, the parameterization (fTU]) is used to obtain 
the IS-RSB free energy as in Eq. ITTI) . However, the pa- 
rameters TUa and TUb are now both taken as variational 
parameters of the free energy 



2F = -JoMaMb + ^ \ + ^^[mp{ql, - q^^) 



f3V^ 
2 



p—a,h 

1 / „ , / „ „^ 1 hi 4 

■p > 



(17) 
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FIG. 6: Order parameters as a function of H for T = 0.0 and FIG. 7: Order parameters versus T for H = 0.4, Jo = Vb = 
two values of T: (a) T = 0.35 and (b) T = 0.3. The disorder 0.95, J = and 1/ = 1.0 with; (a) V = 0.8 and (b) F = 1.0. 
is the same as Fig. ([5]). 



with defined in Eq. p2p . In this approximation to^ ^ 
nib in the asymmetric RSB region. 

The RSB appears in both sublattices simultaneously 
with the same results as those obtained in Eq. pTjl in 
the symmetric region (Ma = —Mb or Ms = 0). How- 
ever, in the asymmetric region {H > and Mg > 0), 
the RSB occurs first at the sublattice b, while the sub- 
lattice a can still present RS solution as it is shown in 
Fig. (O for the particular interaction: Jq — Vq = 0.95, 
J = and V = 1.0 (case (iv)). Nevertheless, both sub- 
lattices present RSB solution at lower temperature when 
r = 0. In Fig. ((5|), the mixed phase is characterized 
when the replica symmetry is broken in the sublattice 
b. However, the spin freezing effects are stronger when 
the two sublattices present RSB solution. For instance, 
the staggered magnetization becomes less dependent on 
T for the region where the sublattices a and b have RSB 
solution (see Fig. Therefore, we can identify two 

distinct asymmetric RSB regions: one with low number 
of frustrated spins where only the set of spins in the sub- 
lattice b exhibits RSB solution; another with the RSB 
solution present in the two sublattices. It occurs because 
the internal fields that act on the sublattices a and b are 
different. In addition, the sublattices can exhibit frustra- 
tion independently (J = 0). 

The free energies given by Eqs. (fTTj) and p7|) are com- 
pared in the insert of Fig. The difference between 
them is very small if it is compared with the difference 
obtained by the RS and the RSB solution with nia = rub 



(see Fig. ([2])). In this sense, we could use the approxima- 
tion m = rua = TTib for the IS-RSB solution with good 
results within the asymmetric RSB region for the case 
with J ^ 0. 

Considering F = 0.0, the i7-dependence of the IS-RSB 
order parameters is exhibited in Fig. ^ for two values 
of temperature: T = 0.35 and T = 0.3. The RSB occurs 
in both sublattices at = 0. When H increases, 5b en- 
hances and 5a decreases in the asymmetric RSB region. 
At r = 0.35 (Fig. (O-a)), 5a becomes zero within an H 
range where only the sublattice 6 presents RSB solution. 
Therefore, the H field can decrease the temperature in 
which the sublattice a exhibits RSB. However, the effects 
of H are not strong enough to produce a stable RS solu- 
tion for a particular sublattice at lower temperatures, as 
it is shown in Fig. ^h) for T = 0.3 where 5a and 5b are 
greater than zero in the mixed phase. 

The effects of F on the order parameters are shown in 
Fig. dZl) for if = 0.4, F = 0.8 and F = 1.0. The RSB re- 
gion of the sublattice a decreases faster than that of the 
sublattice b when F is applied in the asymmetric RSB 
region. Hence, for F sufficiently high (F — 1.0), the RSB 
solution appears only in sublattice 6, while sublattice a 
presents RS solution in all range of T analyzed in Fig. 
^h). Although these results (F > 0) are limited by the 
static approximation, which is expected to yield inaccu- 
rate quantitative results at very low temperature, they 
could provide some qualitative evidences about the RSB 
behavior at least for T > Q}^ For example, when the 
strength of F is enhanced, quantum fluctuations assume 
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a relevant role inducing tunneling between the metastable 
valleys that characterize the SG free energy landscape. It 
could reduce the importance of the RSB, but it is not able 
to restore the RS solution in the sublattice b for instance. 
The quantum fluctuations decrease the region where both 
sublattices present RSB solution. However, they can in- 
crease the region where only sublattice b presents RSB. 
In this case, the quantum fluctuations can decrease the 
number of frustrated sites in the asymmetric RSB region. 

C. Phase Diagrams 

After a detailed study of the asymmetric RSB region 
(mixed phase), we analyze phase diagrams T versus H for 
the two types of coupling previously discussed. The main 
effects of the quantum fluctuations on the phase bound- 
aries and the role of intra-sublattice disordered interac- 
tions on phase diagrams are presented. First, we con- 
sider the thermodynamic behavior of the phase bound- 
aries within a region where the RS solution is stable (at 
T > Tf). In this case, analytical expressions for the free 
energy and the order parameters are achieved by assum- 
ing qp = qpo = Qpi with m = (uip = in Eq. \T^ . The 
transition from the PM to the AF order can be charac- 
terized by the onset of the staggered magnetization Ms . 
If this transition is continuous, we can expand the set 
of saddle-point equations in powers of Ms- After some 
lengthy calculations, we find a compact result: 

Ms = AMs + BMl + CMl + (18) 

where the coefficients A and B are defined in Appendix 
1X1 A second-order phase transition from the PM to the 



AF order is found when the condition A — 1, B <{) and 
C < is obtained. Tricritical points can be found by 
the condition A = 1 and B = with C < 0. We do not 
calculate the coefficient C and so we will not check the 
condition C < 0. 

The phase diagrams for interactions as cases (i), (ii), 
(iu) and (iv) are shown in Figs, (gl-a), dHl-b), dHl-c) and dHl- 
d), respectively. They exhibit a general feature: the AF 
order is destroyed by H at the same time that the mixed 
phase enhances when H increases toward He- Here, He 
is the value of the parallel magnetic field in which Tf 
meets Tjsf. For H > He, the system behaves as in the 
one-lattice model, where Tf decreases monotonically with 
Hfi Particularly, for the case (i), the phase boundaries 
are second-order and the phase diagram (Fig. ([H]-a)) is 
similar to those one obtained in Ref. [l^ for the KS model. 

The role of the intra-sublattice parameters Vq and V 
on the phase diagrams can be analyzed by comparing the 
cases (i) and (ii) (Figs. ([8]-a) and (HJ-b)), and the cases 
(i), (iu) and (iv) (Figs, di-a), di-c) and ©d)), respec- 
tively. In case (ii), Vq is increased at the same time that 
Jo decreases in order to keep r] {Ttq{H = 0)/Tf{H — 0)) 
constant. As a consequence, the AF region reduces and 
a first-order phase transition appears (see Fig dHl-b)). A 
tricritical point is obtained by expansion psp . The first- 
order boundary is located choosing the solution which 
minimizes the free energy. The transition from the RS 
solution to the RSB solution has a discontinuity at He, 
where Tf jumps when it crosses from the region H < He 
to H > He- In this discontinue transition from 

the AF order to the SG phase raises when H increases. 
Particularly for F = 0, the average internal field hp 
{= H + VoMp — JoMp') that acts on the sublattice p is 
directly affected by the relation between Vq and Jq. The 
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hp in which the PM solution becomes stable is achieved 
with smaller H when Vq/Jq enhances. Furthermore, the 
hp for a stable PM solution can be obtained with an H in 
which the AF solution is still energetically favorable (see 
inserts of Fig. (Hl-b)), where the first-order transition 
occurs. There are numerical evidences within IS-RSB 
ansatz which indicate that the mixed/SG phase transi- 
tion can also be first-order (see Fig. (I5]-b)). Therefore, a 
first-order transition can occur when Vq/Jq is sufficiently 
large. 



The intra-sublattice disorder also affects the effective 
AF interaction between the two sublattices as we can see 
by comparing Fig. ^a) with Figs. ([8]-c) and dH-d). For 
instance, the AF region reduces when the V increases at 
the same time that J decreases to keep VV^ + con- 
stant. Thus, a small value of H is able to destroy the 
AF order. Besides, as shown in previous sections, the 
degree of RSB depends on the ration V/J. For example, 
Sb increases with V/ J in the mixed phase when H > 0. 
Consequently, the derivative \dTf/dH\ is higher in the 
AF/mixed phase transition when V increases (see Figs. 
dHl-c) and (Hl-d)). Case (iii) shows a phase diagram simi- 
lar to that one in case (iv). However, the RSB properties 
of these two cases are very distinct as it has been shown 
in sections (|III Ap and (jlll B|) . A first-order phase tran- 
sition also appears for the cases (iii) and (iv). However, 
different from Fig. ([H]-b), the Tf in the asymmetric re- 
gion (AF/mixed phase transition) is higher than in the 
symmetric region (PM/SG transition) for H close to He- 
Consequently, a discontinue transition from the mixed 
phase to the PM phase can occur with the increase of H 
(Figs. ([8]-c) and ([8]-d)). Therefore, the intra-sublattice 
ferromagnetic average interaction can introduce a first- 
order phase transition in this two-sublattice problem. On 
the other hand, the intra-sublattice SK-type of disorder 
changes the behavior of Tf that presents a higher incre- 
ment with H in the AF order when this disorder en- 
hances. 



The internal field is even more complex when F > 0. 
The critical lines and the tricritical points are decreased 
when r increases, as we can see in Fig. ([SJ for phase 
diagrams with F = 1.0. The quantum effects are more 
pronounced at lower temperatures. Particularly, the be- 
havior of Tf is changed by the quantum fluctuations. 
In this case, the derivative \dTf/dH\ is increased by F 
when H is close to He- The region of discontinue tran- 
sition from the AF order to the SG phase (Fig. ^h)) 
and from the mixed phase to the PM phase (Figs. ^ 
c) and (IH]-d)) are also increased by F. This Tf behav- 
ior occurs because of the competition between quantum 
fluctuations (dominant at lower T) and the asymmetry 
caused by H (important near the He) that depends on 
the intra- and inter-sublattices disordered interactions as 
previously shown. 



IV. SUMMARY AND CONCLUSIONS 

We studied a fermionic version for the two-sublattice 
Ising SG model in a transverse (F) and a parallel (H) 
magnetic fields. This model allows SK-type interactions 
for spins of same sublattices and between spins of distinct 
sublattices. The problem is formulated in a Grassmann 
path integral formalism within the static approximation. 
The disorder is treated by the replica trick with Parisi's 
scheme of one-step replica symmetry breaking (IS-RSB). 
A given degree of frustration i] = y/V^ + J^/ {Vq + Jq) 
is used to analyze the main role of each parameter of 
the model. Particularly, we studied the effects of the 
quantum fluctuations on the asymmetric RSB region and 
on phase diagrams. 

For H = 0, the two sublattices show symmetric re- 
sults for the IS-RSB order parameters, which are iden- 
tical for all configuration of intra- and inter-sublattices 
interactions that satisfy the same rj. However, the H 
field introduces an asymmetric RSB region, in which the 
sublattices exhibit different IS-RSB results that depend 
on the configuration of intra- and inter-sublattices dis- 
orders. For instance, the degree of RSB of sublattice b 
(Sb) is higher than Sa- This difference and the value of Sb 
increase when the intra-sublattice disorder {V) enhances. 
Nevertheless, for a general case, the RSB occurs at the 
same temperature for both sublattices when the inter- 
sublattice interactions are disordered (J > 0). In addi- 
tion, the staggered magnetization has a weak dependence 
on T in the mixed phase for IS-RSB solution. Results 
with F > suggest that the quantum tunneling decreases 
the effects of the RSB for both sublattices, but it is not 
able to restore the RS solution at T — > 0. This result is 
in agreement with that one of fermionic one-lattice quan- 
tum SG problem.— 

For a particular case with intra-sublattice disordered 
interactions and only antiferromagnetic coupling between 
the sublattices (J = 0), the RSB appears in each sublat- 
tices independently. The sublattice b can present RSB so- 
lution while the sublattice a still presents RS solution in 
the asymmetric RSB region. However, the ground state 
of both sublattices always occurs with RSB solution for 
F = 0. Nevertheless, the results within IS-RSB suggest 
that the frustration in the mixed phase is changed by the 
F field. In this case, the quantum fluctuations affect with 
more intensity the RSB region of the sublattice a than 
the sublattice 6. The sublattice a can exhibit a stable RS 
solution at any finite temperature with the RSB solution 
present in the sublattice b for a strong enough F within 
the asymmetric RSB region. 

The phase diagrams show a general feature in this 
work: a transition from the PM to the AF order then 
to a RSB region when temperature T decreases for low 
H. The RSB occurs at the freezing temperature Tf with 
a finite staggered magnetization that characterizes 
a mixed phase. A transition from the PM to the SG 
phase (RSB solution with Mg — 0) is observed for high 
H. The H field destroys the AF phase, but it can favor 
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the RSB region within a certain range. In this range, 
Tf increases with H. The intra-sublattice interaction af- 
fects the Tf{H) behavior that presents a higher increase 
with H when the intra-sublattice disorder enhances. The 
presence of an intra-sublattice ferromagnetic long range 
interaction can introduce a first-order phase transition 
from the RS to the RSB solution when H increases. The 
r field always suppresses the magnetic orders. However, 
its effects on Tf{H) are more evident at lower values of 
H, where the asymmetry between the sublattices is small. 



APPENDIX A: EXPANSION FOR THE NEEL 
TEMPERATURE 

The continuous phase transition from the PM phase 
(Ms = 0) to the AF order {Ms > 0) can be obtained by 
expanding the saddle-point equations in powers of Mg. 
After some calculations, we find expansion (|18p where 

A = /3{Vo + Jo)/o'" + f3\V^ - J'){qifS' + ^ifi°) (Al) 
and 

B = P{Vo + Ja)Bi+f3{Vo-Jo)B2+P^{V^-J^)B3, (A2) 
with 

Y-30 

Bi = P^iVo + Jof^ + P{Vo + Jo)P\V^ - J2) 
6 

f21 f20 

x[qi^+xi^]+PHv' - JYki^ifl' (A3) 

fl2 flO 



Xs = [ra - n + qb - qa)/^ = XiM^ + x^M^ H and 

X = {ra+rb + qa + qb)/2 = X2M^ + 0{Mf). In Eqs. 
(lAlll . (|Aa . (lA4l) and ^ 



fF - / Dz 



where i, j, ^ = 0, 1, 2 or 3, and 
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92 



93 



B2 = m2 {f3{Vo + Jo)/o'" + - J') 

x[gi/o" + a;i/^"]}. 



cOi 



S3 = q-3fi^' + x^f^" + P\V' + J')[qi{x2f[ 
+q2C) + xAq2f,' + x^f^)] + P\V' - 

f03 fOO fOl f02 

Mql^+xllj^+q,xl^-+x,qlll-], 



6 



6 



(A4) 



(A5) 



M = {Ma + Mb)/2 ^ m2Af2 + 0(M4), 9, = {qa-qb)/2 = 
qiMs + 93^3 . . . , g = (q, + qb)/2 = q^hP, + 0{Mf), 



K = 



ysinh(Vt/2^^2r2) 



(All) 



with y = (3h and Q defined in Eq. for h = hp in 

the RS solution {Sip = 0) and the PM phase {Ma = Mb, 
qa = qb and = n). 
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